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Abstract. The new data on kt distributions obtained at RHIC are analyzed by means of selected models 
of statistical and stochastic origin in order to estimate their importance in providing new information on 
hadronization process, in particular on the value of the temperature at freeze-out to hadronic phase. 



PACS. 25.75.-q Relativistic heavy ion collisions - 
- 02.50.Ey Stochastic models 



12.40.Ee Statistical (extensive and non-extensive) models 



1 Introduction 

Very recently high kt distributions at RHIC have been re- 
ported in Refs. [1-3]. These data are of potentially high 
interest as a possible source of information on the con- 
ditions existing at the freeze-out in heavy-ion collisions. 
This resulted in a number of works, mostly of statistical 
or thermal origin [4], stressing different possible dynam- 
ical aspects, like the role of resonances or the flow phe- 
nomenon. In our work we would like to show that one can 
account summarily for such (and others) effects consid- 
ered in the literature by using simple minimal extensions 
of the known statistical or stochastic models, which were 
already successfully applied in other analysis of experi- 
mental data. They are: 

(i) The modified statistical model inspired by Tsallis statis- 
tics [5], which generalizes the usual Boltzmann-Gibbs 
statistics to nonextensive systems parametrized by a 
nonextensivity parameter q (for q I one returns to 
the usual Boltzmann-Gibbs extensive scenario); it has 
been already successfully used in this context [6-8] . Pa- 
rameter q summarizes in such approach all deviations 
from the Boltzmann-Gibbs statistics including those 
caused by flow phenomena and resonances [4]. 

{a) A suitable adaptation of the recently proposed model 
derived from the Fokker-Planck equation for the Or- 
stein-Uhlenbeck (0-U) process [9-11] but this time 
used in the transverse rapidity space, i.e., for yt — 
\ ln[(mt + kt)l {mt - kt)] (where = .Jrn?~+Jk^), 
in which one allows for the mass m to be treated as 
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free parameter in order to account for some specific 
features of data (like flow phenomenon) which cannot 
be explained in a usual way. 

As a kind of historical reference point we shall use classi- 
cal statistical model developed long time ago by Hagedorn 
[12] in which transverse momentum distribution of pro- 
duced secondaries is given by the following formula [13] 
(with To being parameter identified with temperature, Th 
denoting the so called Hagedorn temperature [12,13] and 
being pion mass). 



p{m) 



(2) 



As one can see in Fig. 1 although fits to kt distributions at 
^/sWn = 200 GeV obtained by BRAHMS Collaboration 
[1] are quite good, they start to deviate from data at high- 
est values of kt and became very bad there, what is very 
clearly seen in Fig. 2 where we show our fits to STAR data 
[2] covering larger span of transverse momenta. Although 
one can argue that for such large values of kt statistical 
approach must give way to some more detailed dynamical 
calculations [4] , there are examples that suitable modifica- 
tions of statistical approach can lead to quite reasonable 
results in leptonic, hadronic and nuclear collisions. What 
we have in mind here are some non-extensive generaliza- 
tions of statistical model as discussed in [7,6,8]) and some 
specific realization of stochastic approach as proposed by 
[9,10,14]. In what follows we shall therefore apply these 
two methods to nuclear data of Refs. [1-3]. 
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Fig. 1. Results of application of simple statistical model, cf. Eq. (1), to data for fct-distributions at y^ijvjv 
for different centralities by BRAHMS Collaboration [1]. 



200 GeV measured 




k, [GeV/c] 

Fig. 2. Results of application of simple statistical model, cf. 
Eq. (1), to data for fct-distributions at ^sjvjv ~ 200 GeV mea- 
sured for different centralities by STAR Collaboration [2]. 

In next Section we shall analyze data using nonexten- 
sive generalization of statistical model by means of Tsallis 
statistics. In Section 3 we shall analyze data using stochas- 
tic approach in transverse rapidity space. Our conclusions 
are presented in Section 4. 



2 Analysis of kt distributions by generalized 
statistical model based on Tsallis statistics 

In many fields in physics, which use statistical and stochas- 
tic approaches as their tools, it was recognized since some 
time already that the usual Boltzmann-Gibbs approach 
encounters serious problems when applied to systems pos- 
sessing memory effects, correlations (especially long-range 
correlations but also those caused by the production of 
resonances in multiparticle production processes or by the 



Table 1. Values of parameters C, Th and To in eq. (1) used to 
obtain results presented in Figs 1 and 2. The values of x^/n.d.f. 
for BRAHMS data are the same for all centralities and equal to 
19.3/23 and 18.3/23 for mo = 0.5 and 0.55 GeV, respectively. 
For STAR data they are equal 532/32 for C.C = 0-5%, 249/32 
for C.C= 20 - 30% and 308/32 for C.C= 60 - 80%. 



BRAHMS Coll. [1] mo = 0.5 GeV (fixed) 


C.C. 


C 


Th 


To 


(%) 




(GcV) 


(GcV) 


0-10 


177±11 


0.180±0.007 


0.169±0.006 


10-20 


127±9 


0.172±0.008 


0.162±0.006 


20-40 


83±7 


0.156±0.008 


0.149±0.007 


40-60 


44±5 


0.133±0.010 


0.128±0.009 


60-80 


177±11 


0.095±0.0001 


0.093±0.0001 


BRAHMS Coll. [1] mo = 0.55 GoV (fixed) 


0-10 


204±13 


0.172±0.008 


0.162±0.006 


10-20 


146±10 


0.163±0.008 


0.155±0.006 


20-40 


96±8 


0.148±0.008 


0.142±0.007 


40-60 


51±6 


0.124±0.010 


0.121±0.009 


60-80 


204±13 


0.075±0.00007 


0.075±0.0001 




STAR Coll. 


[2] mo = 0.5 GoV (fixed) 


0-5 


816±15 


0.086±0.0001 


0.085±0.0001 


20-30 


382±7 


0.077±0.0001 


0.076±0.0001 


60-80 


106±2 


0.037±0.00001 


0.037±0.00001 



flow effects present there) or which phase space has some 
(multi) fractal structure [5]. Such systems are all, in a 
sense, small, by what we means that the effective range of 
correlations they experience is of the order of dimension 
of the system itself. Therefore they will not show prop- 
erty of extensivity leading to Boltzmann-Gibbs form of 
entropy, which is the basis of any statistical or stochas- 
tical model. One can therefore argue that in such cases 
one has to resort to some dynamical approach in which 
effects mentioned above would be properly accounted for. 
The problem is, however, that there is no unique model 
of this type and usually several approaches are competing 
among themselves in describing experimental data. The 
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other possibility is to realize that most probably our sys- 
tem is not extensive (in the abovementioncd sense) and 
that this fact should be accounted for by using the non- 
extensive form of entropy, for example the so called Tsallis 
[5]. It turns out that such situations are encountered also 
in domain of multiparticle production processes at high 
energy collisions (cf., [6], to which we refer for all details). 
In fact, there already exists a number of detailed analysis 
using a non-extensive approach ranging from kt distribu- 
tions in e+e~ annihilations [7] and in p+p collisions [8] to 
rapidity distributions in some selected reactions [6] . In [7, 
8] a kind of non-cxtcnsivc g- version of Hagcdorn approach 
has been used whereas [6] were exploring information the- 
oretical approach to statistical models including as option 
also its non-extensive version ^ . 

In our work we shall apply Tsallis formalism, treated 
as simplest possible extension of the usual statistical ap- 
proach with parameter q (the so called nonextensivity pa- 
rameter or entropic index) summarizing deviations from 
the usual statistical approach (without, however, specify- 
ing their dynamical origin). It leads to (Tq denotes tem- 
perature): 



2-Kktdkt 



C 



l-(l-9) 



^kf + kf + 
To 



Q 



dki. (3) 



There exist two different formulations leading to slightly 
different forms of parameter Q: 

(a) In first one uses the so-called escort probability distri- 
butions [19], Pi = Pi / ^iPi (cf., for example, analysis 
of kt distributions in e+e" annihilations [7] or in pp 
collisions [8]), in this case Q = q/{l — q). 

(b) In second approach one uses normal definition of prob- 
abilities resulting in Q — 1/(1 — g). In this case, as 
shown in [15,17], parameter q is given by the normal- 
ized variance of all intrinsic fluctuations present in the 
hadronizing system under consideration: 

q = l + u, = l + {{P^)-{Pf)/{f3f. (4) 
This conjecture originates from the observation that: 

-^"'^f{f3)df3 (5) 



[l-{l-q)/3oHo]^ = 

Jo 



where f{0) describes fluctuation of parameter /? and 
has form of Gamma distribution [15,17] (in our case 
Ho = a/ kf + k't + rn? and fluctuations are in tempera- 
ture, i.e., (3 = 1/T and po = (/?) with respect to /(/?)) 



^ It should be mentioned at this point that proper formu- 
lation of Hagedorn model using Tsallis g-statistics has been 
proposed in [16]. We shall not pursue this problem here. 

^ It must be mentioned at this point that this suggestion, 
which in [15] has been derived only for q > 1 case, has been 
shown to be valid also for q < 1 case [17] and extended to 
general form of fluctuations leading then to the new concept of 
superstatistics proposed in [18]. The most recent discussion of 



We have analyzed BRAHMS [1], STAR [2] and PHE- 
NIX [3] data and our results are shown in Fig. 3 and in 
Table 2. It turns out that both form of parameter Q in 
Eq. (3) result in practically identical curves, therefore here 
we are showing only results obtained for Q = g/(l — q). 
The values of parameters are also very close to each other 
with tendency of C, Tq and q estimated by using Q = 
1/(1 — q) being slightly bigger then those obtained for 
Q = q/{l — q). It is worth to stress at this point that such 
comparison of these two approaches has been made for 
the first time here and, as one can see from the presented 
results, it confirms previous expectation (made in [6]) that 
in case of only limited phenomenological applications, as 
is the case of our work, the results from using Eq. (3) with 

Q = q/{l — q) (i.e., parameters: C*^"^ = c, Tq"-* = I and 
qf") — q) are simply connected to those using Q = 1/(1— g) 

(i.e., to parameters: C^''^ = C, r^''^ = L and q^''^ = Q), 
namely: 



Cc^cq. 



(6) 



As one can see from Table 2 these relations are indeed sat- 
isfied (some small difi'erences present could be attributed 
to the fact that both sets of results represent results of sep- 
arate and independent fitting procedures, without making 
use of Eq. (6)). It means therefore that in all phenomeno- 
logical applications one can use either of the two form of 
parameter Q in Eqs. (3), and, if necessary, to use Eq. (6) to 
translate results from one scheme to another. In both cases 
the pion mass value has been used, m = 0.14 GeV (and 
we have checked that additional changes in mass m of the 
type introduced recently in [14], would not affect the final 
results as long as m is limited to, say, m < 0.2 GeV). The 
estimated fluctuations of temperature are of the order of 
30 — 45 MeV. It is interesting to observe that these fluctu- 
ations are weaker at small centralities and grow for more 
peripheral collisions matching very nicely similar fluctu- 
ations seen in p + p data [2] shown here for comparison. 
One should add here also result from similar analysis of 
e+e~ data [7] reporting even higher values of nonexten- 
sivity parameter q (reaching value of g ~ 1-2), i.e., much 
stronger fluctuations. These results confirm therefore, for 
the first time, another expectation made in [6] saying that 
precisely such trend should be observed. This is because 
Eq. (4) can be also interpreted as being a measure of the 
total heat capacity Ch of the hadronizing system (cf. [6]): 



1 

C'h 



1. 



(7) 



As the heat capacity Ch is proportional to the volume, 
Ch ~ V, in our case V would be the volume of interaction 
(or hadronization), it is expected to grow with volume 
and, respectively, q is expected to decrease with V, which 



physical moaning of q parameter when applied to multiparticle 
production processes (and in this context also of the possible 
origin of statistical formulas as well) with relevant references 
can be found in [6]. 
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Fig. 3. Results of application of non-extensive approach given by Eq. (3) with Q = — q) to data for fet-distributions at 
yJSNN = 200 GeV measured for different centralities by BRAHMS [1] and STAR [2] Collaborations and to data at ^/snn = 130 
GeV as measured by PHENIX Collaboration [3]. The results obtained using <5 = 1/(1 — g) instead look essentially the same, 
therefore they are not shown separately. For differences in values of obtained parameters see Table 2. 
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Table 2. Values of characteristic parameters used to fit data on fct-distributious at different ccntrahtics by using non-extensive 
approach as given by Eq. (3) with Q = q/{l — q) and Q = 1/(1 — q) for data at ^/snn = 200 GeV obtained by BRAHMS 
[1] and STAR [2] Collaborations and at ^/sNN = 130 GeV obtained by PHENIX Collaboration [3]. For results obtained using 
Q — 1/(1 — q) we provide also explicit vahies of the corresponding fluctuations of temperature as given by ATo — ^/q — 1 ■ To. 
The order of magnitude of the corresponding errors for To, q and ZiTo: STo, Sq and SATo, , respectively, are listed below as well. 
In analysis wc have used errors either as provided by experiments (for STAR and PHENIX) or assuming systematic error on 
the level of 5% (for BRAHMS). 



BRAHMS Coll. [1] 




Eq. (3) with Q = 
STo =0.005-0.007, Sq 


=0.003-0.005. 


STo 


Eq. (3) with Q = 1/(1 - (?) 
=0.005-0.007, Sq =0.002-0.004, 


















SATo =0.001-0.002 






C.C. 


X7n.d.f. 


C 


To 






xVn-d.f. 


C 


To 


<? 


ATo 


(%) 






(GeV) 










(GeV) 




(GeV) 


0-10 


11.2/23 


1033±78 


0.232 




1.043 


11.2/23 


1033±78 


0.223 


1.041 


0.045 


10-20 


12.9/23 


797±66 


0.224 




1.049 


12.9/23 


797±65 


0.213 


1.047 


0.046 


20-40 


12.7/23 


525±49 


0.215 




1.055 


12.7/23 


525±49 


0.204 


1.052 


0.047 


40-60 


10.5/23 


304±38 


0.193 




1.067 


10.5/23 


304±38 


0.181 


1.063 


0.045 


60-80 


2.85/22 


41±5 


0.175 




1.084 


2.85/22 


41±5 


0.161 


1.077 


0.045 


STAR Coll. [2] 




Eq. (3) with Q = 










Eq. (3) with 


Q = V(i 


-q) 






STo =0.002-0.003, Sq 


= 0.001-0.002. 


STo =0.002-0.003, Sq ^0.001, 


















SATo 


5^0.001. 






C.C. 


xVn.d.f. 


C 


To 




g 


xVn.d.f. 


C 


To 


q 


ATo 


(%) 






(GeV) 










(GeV) 




(GeV) 


0-5 


170/32 


4684±231 


0.171 




1.071 


170/32 


4686±231 


0.159 


1.066 


0.041 


5-10 


68/32 


3393±184 


0.176 




1.068 


67.8/32 


3393±185 


0.165 


1.064 


0.041 


10-20 


69/32 


2767±144 


0.171 




1.073 


69.2/32 


2768±144 


0.160 


1.068 


0.042 


20-30 


45/32 


1928±102 


0.169 




1.075 


44.7/32 


1928±102 


0.157 


1.070 


0.042 


30-40 


44/32 


1391±78 


0.165 




1.078 


43.9/32 


1391±78 


0.153 


1.072 


0.041 


40-60 


19/32 


896±50 


0.153 




1.085 


19.2/32 


896±50 


0.141 


1.079 


0.040 


60-80 


14/32 


414±25 


0.137 




1.095 


14.2/32 


413±25 


0.125 


1.087 


0.037 


1> + P 


9.7/29 


02=7 


().il7 




1.099 


l).()2/29 


(ii. 9=7.1 


0.107 


1.090 


().0.:i2 


PHENIX Coll. [3] 




Eq. (3) with Q = 
STo =0.011-0.016, Sq 


a/(i-g) 

=0.008-0.010. 


STo 


Eq. (3) with Q = 1/(1 -q) 
=0.011-0.016, Sq =0.005-0.011, 


















SATo =0.003-0.005 






C.C. 


xVn.d.f. 


C 


To 




g 


xVn.d.f. 


C 


To 


q 


ATo 


(%) 






(GeV) 










(GeV) 




(GeV) 


0-5 


5.13/29 


1694±409 


0.201 




1.049 


5.1/29 


1694±411 


0.192 


1.047 


0.042 


5-15 


3.62/29 


1330±316 


0.199 




1.051 


3.6/29 


1330±316 


0.190 


1.048 


0.042 


15-30 


5.55/29 


846±206 


0.196 




1.054 


5.6/29 


846±206 


0.186 


1.051 


0.042 


30-60 


2.63/29 


433±113 


0.178 




1.066 


2.6/29 


433±113 


0.167 


1.074 


0.045 


60-80 


10.6/29 


139±48 


0.152 




1.080 


10.6/29 


139±48 


0.141 


1.062 


0.035 


80-92 


9.10/29 


74±45 


0.121 




1.098 


9.1/29 


74±42 


0.110 


1.089 


0.033 



is indeed the case if one puts together results for e+e ,pp 
and A A coUisions. 



3 Analysis of kt distributions using stochastic 
approach in yt space 

Whereas previous approach was concerned with extension 
of the purely statistical approach the one presented now 
will go bit further by modelling hadronization process by 
a kind of diffusion mechanism Refs. [9,10,14] in which 
the original energy of projectiles is being dissipated in 
some well defined way into a number of produced sec- 



ondaries occurring in different part of the phase space ^. 
In the case considered here it is diffusion process taking 
place in the kt space. Actually, it turns out that it is more 
suitable to consider such diffusion as taking place in the 
yt = smh~^ (kt/m) space. In this case one obtains the fol- 
lowing Fokker-Planck equation: 



dPt{yt ,t) 
dt 



= 7 



dytPt{yt, t) ^a^d^Pt{yt,t) 



dyt 



27 



dy! 



(8) 



It should be mentioned here that there exist also non- 
extensive versions of such diffusion process applied to mul- 
tiparticle production data [20] but we shall not pursue this 
possibility here. 
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Table 3. Values of characteristic parameters used to fit data 
on fet-distributions at different centralities by using stochastic 
approach as given by Eq. (9) and presented in Fig. 4 for data 

at yi]^ = 200 GeV obtained by BRAHMS [1] and STAR [2] 
Collaborations and at ^/snn = 130 GcV obtained by PliENIX 
Collaboration [3] . The order of magnitude of the corresponding 
errors for To, 5To, are listed below as well. 



BRAHMS Coll.; 5Tb = 0.008-0.012; Sm = 0.024-0.031 



C. C. (%) 


xVn.d.f. 


C 


To (GeV) 


m (GeV) 


0-10 


39.9/23 


140±9 


0.201 


0.784 


10-20 


24.2/23 


108±7 


0.199 


0.725 


20-40 


20.1/23 


72±5 


0.196 


0.671 


40-60 


11.9/23 


38±4 


0.185 


0.577 


60-80 


4.06/22 


4.3±0.5 


0.184 


0.515 


STAR Coll.; STo = 


0.004-0.006; 5m = 0.009-0.014 


C. C. (%) 


xVn.d.f. 


C 


To (GeV) 


m (GeV) 


0-5 


221/32 


484±22 


0.169 


0.533 


5-10 


124/32 


370±18 


0.170 


0.547 


10-20 


121/32 


310±14 


0.168 


0.513 


20-30 


92.9/32 


217±10 


0.168 


0.498 


30-40 


89.9/32 


158±8 


0.165 


0.473 


40-60 


43.9/32 


99.6±5.0 


0.157 


0.419 


60-80 


22.3/32 


43.2±2.3 


0.143 


0.349 


P + P 


17.4/29 


5.29±0.63 


0.126 


0.298 


PHENIX Coll.; STo 


= 0.020-0.037; 5m = 0.058-0.078 


C. C. (%) 


xVn.d.f. 


C 


To (GeV) 


m (GeV) 


0-5 


8.06/29 


161±34 


0.185 


0.734 


5-15 


5.61/29 


124±26 


0.185 


0.729 


15-30 


7.27/29 


80±18 


0.185 


0.700 


30-60 


3.49/29 


39±9 


0.177 


0.594 


60-80 


11.1/29 


12±4 


0.158 


0.460 


80-92 


9.01/29 


6.1±3.6 


0.131 


0.317 



Its solution can be expressed by Gaussian distribution, 



2'Kktdkf. 



CPt{yt,t) 



C 



exp 



with 



1 



-27* 



(9) 



(10) 



In Fig. 4 wc show our results of using Eq. (9). It should 
be noticed that now, following [14] , we have regarded mass 
m to be a free parameter. Only then we can obtain good 
agreement with data. In a sense, variable mass m corre- 
sponds in this approach to the non-extensivity parameter 
q introduced in Section 2 in that it summarily accounts for 
some additional effects not accounted by simple diffusion 
process (like, for example, effect of resonances and flow). 

In the stochastic approach considered here we do not 
have direct access to the temperature Tq. It is accessible 
only if we additionally assume validity of the Einstein's 



See also [9]. Actually, Eq. (9) is the same as the formula 
used already long time ago in [11]. 



fluctuation-dissipation relation, which in our case means 
that measure of the size of diffusion (dissipation), V^{t), 
can be expressed by the temperature To and mass m: 



m 



(11) 



Therefore our results for V'^ shown in Fig. 4 (see inlets), 
where V{t)'^ increases with increasing centrality, would 
indicate that temperature Tq obtained by applying Ein- 
stein's relation with m kept constant would increases with 
centrality as well, contrary to what has been obtained 
above by applying q-statistical approach. We have allowed 
then (following [14]) the mass m to be a free parame- 
ter and the best fit is obtained when m decreases with 
centrality, see inlets in Fig. 4. The resulting temperature. 
To ~ m • V(2, behaves then in essentially the same way as 
function of centrality as in the g-statistical approach, cf.. 
Table 3 and Fig. 5^. 



4 Concluding remarks 

We have provided here systematic analysis of recent RHIC 
data on kt distributions [1-3] by using three different kinds 
of statistical approaches: Hagedorn model [12], two ver- 
sions of the modified statistical based on Tsallis statistics 
[5] and a suitable adaptation of the stochastic model pro- 
posed in [9]. We have found that Hagedorn-type model 
[12] cannot fit data at large kt (its widely used for quick 
estimations simplified version with p{m) = 1, which is 
then just a simple Boltzmann gas model with only one 
parameter, the temperature Tq, fails completely even for 
smaller kt, cf. Table 4). However, these data can still be 
reasonably well fitted either by non-extensive extensions 
of statistical model [5] or by picture of some suitable dif- 
fusion process taking place in transverse rapidity space 
[9,10]. This is specially true if one limits itself to fc^ < 5 
GeV/c range as the case of BRAHMS [1] and PHENIX [3] 
data, the kt < 12 GeV/c range considered in STAR exper- 
iment [2] seems to be already too big to be fitted properly 
even with these two approaches (the corresponding values 
of the Suie considerably bigger in this case and the val- 
ues of parameters obtained for STAR and BRAHMS data, 
which were taken at the same collision and at the same 
energy, are also different). 

As is shown in Fig. 5, the temperatures To obtained in 
modified statistical and stochastic approaches (with vary- 
ing mass m) are essentially very similar to each other and 
follow the same dependence on the centrality, namely Tq 



* It should be stressed hero that for constant value of mass, 
m = 0.14 GeV as used for g-statistics case above, we would 
have obtained somewhat higher values of x^'s- In addition, it 
is interesting to observe at this point that the fact that we can 
fit data within modified stochastic approach only by allowing 
for a kind of "quasiparticlcs" of mass m, different for different 
centralities, corresponds in a sense to introducing parameter q 
to the usual statistical model. The possible dynamical origin 
and meaning of such variable mass is, however, still lacking. 
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Fig. 4. Results of application of stochastic approach as given by Eq. (9) to data for fct-distributions at ^ysNN = 200 GoV 
measured for different centralities by BRAHMS [1] and STAR [2] Collaborations and at ^/s]^ = 130 GeV obtained by PHENIX 
Collaboration [3]. Notice that mass m is treated here as free parameter, in similar way as in [14]. 
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Table 4. Comparison of investigated models: simple statistical model (i.e., Hagedorn model as given by eq. (1) but with 
p(m) = 1, in which case it is just a simple statistical Boltzmann gas model with only one parameter, namely temperature To), 
non-extensive Tsallis distribution (NETD) and Ornstein-Uhlenbeck process (0-U), using data on kt distributions at s/snn = 200 
GeV obtained by BRAHMS Collaboration [1] for smallest and largest centralities. 



C.C (%) 


Simple statistical model, 
Eq. (1) with p{m) — 1 
To q rn 

(GeV) (GeV) 


NETD Eq. (3) 
(with g = g/(l - q)) 
To q m 
(GeV) (GeV) 


O-U Eq. (9) 

To q m 
(GeV) (GeV) 


0-10 X^'/n.d.f 


177/23 
0.302 — — 


10.2/23 
0.232 1.043 — 


39.9/23 
0.201 — 0.784 


60-80 xVn.d.f 


567/22 
0.325 — — 


2.76/22 
0.175 1.084 — 


4.06/22 
0.184 — 0.515 
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Fig. 5. Comparison of temperatures of hadronization obtained 
by using different approaches as given by: (a) - Eq. (3) with 
Q = q/{l - q); (6) - Eq. (3) and Q = 1/(1 - q); (c) - Eq. 
(9). In the later case To has been obtained from the values 
of V^^ and m obtained in Fig. 4 by using Einstein's relation: 
To =m- ^^(t). 



the two possible approaches using 5-statistics are equiv- 
alent to each other, at least in the frame of Umitcd phe- 
nomenological approach presented here. One should also 
stress at this point that q statistical approach offers unique 
information on fluctuations in the system, which can be 
translated into information on its volume. Our results for 
AA and pp collisions taken together with old results for 
e+e" annihilations indicate in this respect distinct growth 
of the expected volume of interactions from the most ele- 
mentary annihilation processes to the nuclear collisions. 
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decreases when collision is more peripheral. However, be- 
cause stochastic approach seems to be more dynamical 
than g-statistical one (where the true dynamical origin of 
the nonextensivity parameter is not yet firmly established, 
see [6, 15]), we regard as the most valuable our finding that 
stochastic approach [9,10,14] works so well and can serve 
to provide first simple estimations of any new data in the 
future. On the other hand we have also demonstrated that 
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